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Objective 

• Validate the newly proposed nonlinear turbulence models for 
momentum and scalar transport 

• Evaluate the newly proposed scalar joint probability density functions 
(APDF and DWFDF) along with its Eulerian method in the National 
Combustion Code (NCC). 

• Simulations conducted include 

- Steady Reynolds averaged Navier-Stokes RANS, 

- Unsteady RANS (URANS) 

- Time-filtered Navier-Stokes (TFNS) --- very large eddy simulation 

- Hybrid RANS/APDF 

- Hybrid URAND/APDF 

- Hybrid TFNS/DWFDF --- very large eddy simulation 

In the hybrid scheme, the transport equations of energy and species are 
replaced by the APDF or DWFDF equation 

Some positive effects of nonlinear models and hybrid 
approaches observed. 
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Confined Swirling Coaxial Jets 


Geometry configuration 
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Computational domain and grid 

849,189 tetrahedral elements 
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Simulation of water jet using NCC 


• Experiments are water jets. 

• NCC code is for ideal gas flow. 

• “Reynolds number similarity law” under low speed 
conditions was used for rescaling between the water 
and gas flows. 


3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
URANS-s p_L, 2nd = -.01 , 4th=.05 
dt=1.e-05, 18,000 time step. 
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Outline 

• Basic equations for RANS, URANS and 
TFNS. 

• Scalar APDF and DWFDF equation for 
hybrid approach. 

• Comparison of numerical simulations 
with experimental data. 

• Conclusions. 
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Equations for RANS, URANS and TFNS 
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Nonlinear models 


p UU . , pU i e , p U i O m can be expressed through the following 


turbulent transports (stresses and scalar fluxes) T^- , © ? : 




* J 


i J 


Nonlinear models are (NASA/TM-1997-113112, 2010-216323): 

r v =\s v t a -2C M -fp^{S 9 -S v Sj3)-A i -f-p^{s^-O ik S tl ) 
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where f - 0 < f(RCP) < 1, 0 < PCP < 1 Resolution control parameter. 
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Equation for Scalar APDF & DWFDF, ; x,t) 

A density weighted ensample averaged or time filtered 
fine grained scalar PDF, F n> (y/ ; x,t), is defined as: 

F^iy/ \ x ,t) = p(x,t')f^(y/;x,t')G(t-t')dt' 

where, (ys ; x , t ') = 5 ( O ( x, t ’ ) - y / ) — fine grained scalar PDF 

With above definition, the density weighted mean or filtered variable O 
can be fully expressed using (y / ; x, t) as: 

[ i/s F^iys ; x,t)dys = p<l> = p(x,t)<S>(x,t) 

Insert this relationship into scaler transport equation, we will 
obtain a transport equation for scalar F (l> (y / ; x,t ) : 
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The resulting transport equation for F a> (y/:x,t) is 
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Note that the chemical reaction term S k (y) is in a closed form! 


• The diffusion term in the sample space is further simplified 
to fit the available PDF solution procedure built in NCC 
code as 
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Results of simulations and comparisons 
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Global flow features of TFNS (VLES) simulation 


Scalar flux model: linear 


Scalar flux model: non-linear 


3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
TFNS with sp_L, wf, 2nd=0, 4th=.05 
dt=2.e-06, 335,000 time step. 
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3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
TFNS with sp_nl, wf, 2nd=0, 4th=.05 
dt=2.e-06, 335,000 time step. 
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Time history of velocity components u, v, w at 

probe 4 


Scalar flux model: linear 


TFNS-sp_L 



Scalar flux model: non-linear 


TFNS-sp_nl 
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Vortex brake-down Bubble 


Scalar flux model: linear Scalar flux model: non-linear 


3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
TFNS with sp_L, wf, 2nd=0, 4th=.05 
dt=2.e-06, at 335,000 time step 
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3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
TFNS with sp_nl, wf, 2nd=0, 4th=.05 
dt=2.e-06, at 335,000 time step 
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Iso-surface of U=0 




Iso-surface of U=0 
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Contour of instant axial velocity, u 


Scalar flux model: linear 


Scalar flux model: non-linear 


3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
TFNS with sp_L, wf, 2nd=0, 4th=.05 
dt=2.e-06, at 335,000 time step 
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3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
TFNS with sp_nl, wf, 2nd=0, 4th=.05 
dt=2.e-06, at 335,000 time step 
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Contour of colored 02 concentration 


Scalar flux model: linear Scalar flux model: non-linear 


3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
TFNS with sp_L, wf, 2nd=0, 4th=.05 
dt=2.e-06, at 335,000 time step 
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3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
TFNS with sp_nl, wf, 2nd=0, 4th=.05 
dt=2.e-06, at 335,000 time step 
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Contour of effective viscosity, 


mu 


Scalar flux model: linear 


Scalar flux model: non-linear 


3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
TFNS with sp_L, wf, 2nd=0, 4th=.05 
dt=2.e-06, at 335,000 time step 
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3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
TFNS with sp_nl, wf, 2nd=0, 4th=.05 
dt=2.e-06, at 335,000 time step 
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Contour of vorticity magnitude 


Scalar flux model: linear 


Scalar flux model: non-linear 


3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
TFNS with sp_L, wf, 2nd=0, 4th=.05 
dt=2.e-06, at 335,000 time step y 

a 


Vorticity Magnitude 

3 .000E+03 
2 .667E+03 
2 . 333E+03 
2 .000E+03 
1.667E+03 
1 „ 333E+03 
1.000E+03 
6 .667E+02 
3 . 333E+0 2 
0 .000E+00 



3D Roback-Johnson flow 
Inner-dye-jet, Turb Intens = 0.01 
TFNS with sp_nl, wf, 2nd=0, 4th=.05 
dt=2.e-06, at 335,000 time step y 
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U (m/s) 


URANS simulation 




URANS: Axial velocity along centerline Inner Jet concentration at downstream 
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U (m/s) 


TFNS (VLES) simulation 




TFNS: Axial velocity along centerline Inner Jet concentration at downstream 
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Inner Jet Concentration 


Hybrid TFNS/DWFDF vs. TFNS 


TFNS, dt=2.e-06, 2nd=0., 4th=0.05, cfl=1.0 



TFNS, dt=2.e-06, 2nd=0., 4th=0.05, cfl=1.0 



Mean concentration and axial velocity distribution along the centerline 
Positive improvements shown from hybrid method 
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Appendix: RANS simulation 
standard k — s model vs. nonlinear model 


RANS standard k-eps model with WF 
cfi=1.0, 2nd=-.01 , 4th=.Q5, conv=1.e-03 
converged at 155,224 iteration 
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3D Roback-John son flow 
Inner-dye-jet, Turb intens = 0.01 
RANS-NL-wf, cfi=1.0, 2nd=-.01 , 4th= 05 
At 153,224 iteration, converged Y 
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Center recirculation zone missed 


Center recirculation zone 
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Appendix: RANS simulation 
standard k — s model vs. nonlinear model 


RANS standard k-eps model with WF 
cfl=1 .0, 2nd=-.Q1 , 4th=.G5, conv=1.e-03 
converged at 155,224 iteration 



x 


02red 

mm 2 . 3 DOE -01 
I 2 . 0 44E-01 
1.789E-01 
1 * 533E-01 
1.278E-01 
1.G22E-01 
7 * 6 67E-02 
5 . 11 IE -02 
J 2.556E-02 
■ O.OOOE+OO 


3D Roback-Johnson flow 


Inner-dye-jet, Turb Intens = 0.01 
RANS-NL-wf, cfl=1.0, 2nd=-.01 ( 4th-. 05 
At 153,224 iteration, converged 
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Contour of inner jet concentration at center plane (x,y) 
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U (m/s) 


Appendix: RANS simulation 
standard k — s model vs. nonlinear model 


RANS, R-J coaxial jet flow 




Axial velocity and concentration along the centerline 
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Conclusions 

• Two groups of validations have been performed against experimental 
data: 

- The first group focuses on the turbulent scalar flux models: linear vs. 
nonlinear. Simulations include RANS, URANS and TFNS. 

- The second group focuses on the hybrid approach. Simulations include 
RANS/APDF, URANS/APDF and TFNS/DWFDF. 

• Regarding the scalar flux model: 

- the linear and nonlinear scalar flux models have the same or similar 
behaves in RANS, URANS and TFNS simulations. 

- In the case of TFNS simulation, TFNS results demonstrate significant 
improvements over their RANS and URANS counter parts. 

• Regarding the hybrid approach: 

- RANS/APDF, URANS/APDF and TFNS/DWFDF simulations show that they 
are quite close to their respective RANS, URANS and TFNS counterpart. 

- The hybrid approach appears to be more robust in the unsteady 
calculations and converge faster to use less computing time. 

- The above observations show a quite positive opinion of present hybrid 
PDF method for even non-reacting flow simulations. 
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